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In hole-doped high-T c superconductors, thermal conductivity ft increases drastically just below T c , 
which has been considered as a hallmark of a nodal gap. In contrast, such a coherence peak in « is 
not visible in electron-doped compounds, which may indicate a full-gap state such as a d + is-wave 
state. To settle this problem, we study k in the Hubbard model using the fluctuation-exchange 
(FLEX) approximation, which predicts that the nodal d-wave state is realized in both hole-doped 
and electron-doped compounds. The contrasting behavior of k in both compounds originates from 
the differences in the hot/cold spot structure. In general, a prominent coherence peak in k appears 
in line-node superconductors only when the cold spot exists on the nodal line. 
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mation 



In strongly correlated electron systems, transport phe- 
nomena give us significant information on the many- 
body electronic states. In high-T c superconductors 
(HTSCs), for example, both the Hall coefficient i?H 
and the thermoelectric power S are positive in hole- 
doped compounds such as YBa 2 Cu307_a (YBCO) and 
La2-5Sr5Cu04 (LSCO), whereas they are negative in 
electron-doped compounds like Nd2-<5Ce,sCu04 (NCCO) 
and Pr2-<5Ce,5Cu04 (PCCO) [1]. These experimental 
facts originate from the difference in the "cold-spot," 
which is the portion of the Fermi surface where the re- 
laxation time of a quasiparticlc (QP), t^, takes the max- 
imum value [2-4] : 

Below T c , electronic thermal conductivity k has been 
observed intensively since it gives us considerable infor- 
mation on the superconducting state; it is the only trans- 
port coefficient which remains finite below T c . For exam- 
ple, the k-dependence of the SC gap can be determined 
by the angle resolved measurement of k under the mag- 
netic field [5, 6]. Also, one can detect the type of nodal 
gap structure (full-gap, line-node, or point-node) by mea- 
suring k at low temperatures (T <C T c ). For T < T c , k 
also shows rich variety of behavior in various supercon- 
ductors. In conventional full-gap s-wave superconduc- 
tors, the opening of the SC gap rapidly decreases the den- 
sity of thermally excited QPs, causing k to decrease. On 
the other hand, k shows "coherence peak" behavior just 
below T c in several unconventional superconductors with 
line-node gaps, e.g., hole-doped HTSC [7-9], CeCoIn 5 
[10, 11], and URu 2 Si2 [12]. A previous theoretical study 
based on a BCS model with c?-wave pairing interaction 
[13] discussed that the coherence peak in YBCO origi- 
nates from the steep reduction in r below T c . 

In sharp contrast, no coherence peak in k is observed 
in electron-doped HTSCs [14, 15], irrespective that a re- 
cent ARPES measurement [16] suggests that the d x 2. y 2- 
wave state is realized. The observed k-dependence of 
the SC gap function in NCCO, which prominently de- 
viates from cos fez — cosfcj,, is well reproduced by the 



fluctuation-exchange (FLEX) approximation [17], which 
is a self-consistent spin fluctuation theory. On the other 
hand, recent point-contact spectroscopy for PCCO [18] 
suggests that a full-gap SC state such as d x 2_ y 2 + is or 
d x 2_ y 2 + id xy state is realized for 6 — 0.15 and 0.17. To 
find out the real SC state in electron-doped HTSC, we 
have to elucidate whether the "absence of coherence peak 
in k" is a crucial hallmark of the full-gap SC state, or it 
can occur even in nodal gap superconductors. 

In this letter, we present a theoretical study of the elec- 
tronic thermal conductivity k in HTSCs using the FLEX 
approximation. This is the first numerical study of trans- 
port properties in the SC state based on the repulsive 
Hubbard model. In deriving the relaxation time Tk, both 
the strong inelastic scattering due to Coulomb interaction 
and weak elastic impurity scattering are taken into con- 
sideration, which corresponds to optimally-doped YBCO 
and NCCO samples, respectively. We find that a sizable 
coherence peak of k in YBCO originates from the reduc- 
tion in inelastic scattering. In contrast, the coherence 
peak is absent in NCCO in spite of that d x 2_ y 2-w&ve SC 
state is realized, since the nodal point does not coincide 
with the cold spot in the normal state. Thus, contrast- 
ing behaviors of k in YBCO and NCCO are explained on 
the same footing as d x 2_ y 2-w&ve superconductors. This 
result was not derived in the BCS model [13]. 

Here, we study the following repulsive Hubbard model: 

H = H £ k c L c ka + jj c k+qT c k'-qi c k'|C kT (1) 

k.cr k,k'q 

where U is the Coulomb interaction and 
£k = 2t (cos(k x ) + cos(k y j) + 4ti cos(k x ) cos(k y ) + 
2t2 (cos^A^ ) + cos(2fcj,)) is the kinetic energy of free 
electrons. Hereafter, we put t = —1.0, t\ = 0.167, 
and ti = —0.2 [2] to reproduce the Fermi surface of 
YBCO and NCCO. We also put U = 8.0 for YBCO and 
U = 5.4 for NCCO. Here, no phenomcnological fitting 
parameters are introduced except for U. 
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FIG. 1: Temperature dependence of the DOS given by the 
FLEX approximation for YBCO (n = 0.9). Note that T c = 
0.024. 



FIG. 2: Position of the hot/cold spots for YBCO (n = 0.9) 
and for NCCO (n = 1.10). AFBZ represents the antiferro- 
magnetic Brillouin zone. 



In the FLEX approximation, the normal and anoma- 
lous self-energies are given by 

E n (k,ie n ) = — ^r- }] Gg +k (ie n + icjj) 
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where e n = 7rT(2n + 1) and uji — 2nTl are the Matsub- 
ara frequencies for fermions and bosons, respectively. Xs 
and Xc are the dynamical spin and charge susceptibilities, 
which are given by 
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where G and i* 1 are the normal and anomalous Green 
function, respectively. They are given by 

G{ie n ) = (i e „ + f : k + S"(- J e„)) J D(e n )- 1 (8) 
F{ie n ) = Y> a {-ie n )D(ie n )~ l (9) 
D(ie n ) = (-e n -e k + £ n H£„))(-e„ + g k + £"(ie„)) 
-(E a (-ie n )) 2 (10) 

where e k = e k — fi; /j is the chemical potential. In 
the FLEX approximation, we solve eqs. (2)-(10) self- 
consistently by choosing /i to adjust the electron filling n. 



In the following numerical study, we use 64 x 64 k-meshes 
and 2048 Matsubara frequencies. Figure 1 represents the 
density of states (DOS); p(ui) — jf J2 k i^) / 71 ■ Here, 
the advanced (retarded) Green function G k {ui) (G k (u>)) 
is given by the numerical analytic continuation of the 
Matsubara Green function from the lower (upper) half 
plane in the complex lu space. 

Figure 2 shows the location of the hot/cold spots for 
both YBCO and NCCO in the normal state. The trans- 
port phenomena are governed by QPs around the cold 
spot, where the QP damping rate 7 k = Im£ k (— i8) 
[= l/2r k ] takes the minimum value. According to 
the FLEX approximation, the cold-spot in hole-doped 
[electron-doped] systems is around (7r/2,7r/2) [(n, 0)] [2]. 
The position of the cold-spot in electron-doped systems 
was confirmed by ARPES measurements [19, 20] after 
the theoretical prediction [2] . 

Hereafter, we derive the electric thermal conductivity 
k. According to the linear response theory [21-23], 

« = ^ f dr V l(Q fci (r)Q fe2 (0))e— (11) 



Qkx{r) = q kx ] + q kxl + q£ x + q£ x (12) 

where Q kx is the heat current operator in the supercon- 
ducting state: q kxa and q^ x are given by 



qkxa = uJv kxa , q£ x = LUV kx , 



(13) 



where v kxa represents the "Fermi velocity" and vfc x is the 
"gap velocity" [21]: 

Vkxa = v kx cl a c kt7 , v£ x = uLc k T c -k|> (14) 

n dY? , <9E? 
where v kx = v Ux + — and v£ x = — . 

As a result, the expression for k in the SC state with 
dropping the current vertex correction (CVC) is given by 
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FIG. 3: Schematic T-dependences of 7k LEX and 7i mp . 
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where f(z) = (e z / T + 1) — 1 . In the normal state, heat 
CVC due to Coulomb interaction is small, as shown in the 
second-order perturbation theory with respect to U [23] , 
and in the FLEX approximation [25]. This fact will also 
be true below T c since the particle-nonconserving four 
point vertex is much smaller than the particle-conserving 
one [26]. Therefore, we neglect the heat CVC in the 
present numerical study. We find that the first term in 
eq. (15), which is proportional to {qk x } 2 , is predomi- 
nant, and the other terms which contain the gap velocity, 
g£ x , are negligibly small. On the other hand, the charge 
CVC gives anomalous transport properties for i?n [2], 
S [3], magnetoresistance [24] and Nernst coefficient [25] 
both in hole-doped and electron-doped systems. Anoma- 
lous transport phenomena due to charge CVC are also 
observed in CeMIn 5 (M=Co,Rh) [27] and in re-(BEDT- 
TTF) 2 X [28]. 

Here, we include the QP damping due to impurity scat- 
tering 7; mp by replacing E raii (k, u) — > S nit (k, w) 
Then, the total QP damping rate is 7k 
[13, 29], where 7£ LEX = ImS nR {k,uj). In the i-matrix 
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where nj mp is the impurity concentration, 
purity potential, and go = -h ^ k is the local Green 
function. Hereafter, we consider the nearly unitary limit 
case I ~ 00 where 7i m p(T = 0) takes a constant value in 
the self-consistent calculation [13, 29], and assume that 
7imp < 7k at T = T c . In this case, we are allowed 
to put 7imp(T) = 7imp(T = 0) since the T-dependence 
of 7imp affects K near T c only slightly. A schematic In- 
dependences of 7k LEX and 7i mp are shown in Fig. 3. 

Figure 4 represents the temperature-dependence of k 
given by eq. (15). In YBCO, k increases drastically 
below T c since the AF fluctuations, which are the origin 
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FIG. 4: Obtained k for YBCO (n = 0.9; T c = 0.024) and 
for NCCO (n = 1.15; T c = 0.014), for 7imp = 0, 0.002 and 
0.004. "Normal" represents k given by the normal state FLEX 
approximation under the constraint F = E a = 0. 



of inelastic scattering, are reduced due to the SC gap. 
Since 7k is much larger than 7i mp at T > T Cl k in the 
normal state is affected by 7i mp only slightly. For T <C T c , 
on the other hand, k is suppressed by 7i m p; k shows the 
maximum when 7k LEX ~ 7imp is satisfied at the nodal 
point. The obtained result is consistent with experiments 
[7-9]. In strong contrast, in NCCO, "coherence peak" in 
k is very small even for 7i mp = 0, which is also consistent 
with experiments [14, 15]. 

Here, we discuss the reason why the coherence peak 
in k is present in YBCO whereas it is absent in NCCO. 
Below T c , only thermally excited QPs above the SC gap 
can contribute to re, except at the nodal point. According 
to eq. (15), thermal conductivities in the normal state 
(re„) and in the line-node SC state (re s ), where p(e) oc |e|, 
are approximately given by 



re n oc 

K s OC T 2 /7nodc 



(16) 
(17) 
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FIG. 5: T-dependence of 7 k for both YBCO and NCCO for 
7im P = 0. Thick full lines and thick broken lines represent 
Tnodo (T < T c ) and 7 co i d (T > T c ), respectively. 




YBCO 



T 

FIG. 6: n-dependence of k for YBCO; n = 0.8 ~ 0.9. 

where 7 co id and 7 no de represent 7k at the cold spot above 
T c and that at the nodal point below T c , respectively. 

Figure 5 shows the T-dependencc of 7k given by the 
FLEX approximation. In YBCO, both 7 co id and 7 no do 
are given by 7k at the same point; k w (tt/2, tt/2). As the 
temperature drops, 7 co id decreases moderately in propor- 
tion to T in the normal state. Below T c , 7 no( je quickly 
approaches zero since inelastic scattering is suppressed 
by the SC gap. As a result, n shows a prominent co- 
herence peak below T c , as recognized by eqs. (16) and 
(17). In NCCO, on the other hand, 7 co id is given by 7k 
at k « (it, 0), which is different from the nodal point of 
the SC gap; k w (tt/2, tt/2). Since 7 no do is much larger 
than 7 co id at T = T c in NCCO as shown in Fig. 5, k is 
not enhanced in NCCO below T c . Although the numeri- 
cal accuracy becomes worse for NCCO below T ~ 0.015, 
the obtained result of k will be qualitatively reliable. In 
summary, the coherence peak in k is absent even in nodal 
SC when the cold spot and the nodal point are different. 

Figure 6 shows the obtained doping dependence of k in 



YBCO. In over-doped case (n = 0.8), the enhancement 
of k is largest, and it decreases in optimally (n = 0.85) 
and under-doped (n = 0.9) cases since the cold spot ap- 
proaches the AFBZ as n — > 1. This tendency is consistent 
with experiments [7]. Note that in real materials, T c in 
under-doped case is smaller than that in optically-doped 
case. In the FLEX approximation, however, T c mono- 
tonically increases as n — > 1 since the pseudo-gap state 
in under-doped region cannot be described. The char- 
acteristic pseudo-gap phenomena are well described by 
including the self-energy correction due to strong SC fluc- 
tuations into the FLEX approximation, which is called 
the FLEX+T- matrix approximation [25, 30]. 

In the present work, we assumed that the inelastic scat- 
tering is dominant, and neglected the temperature depen- 
dence of 7imp- This assumption will be allowed for clean 
optimally-doped HTSCs. In dirty samples where elastic 
scattering is large, we should calculate the T-dependcncc 
°f 7imp using the self-consistent i-matrix approximation 
[13, 29]. In under-doped systems, however, the t-matrix 
approximation is not sufficient since the radius of "ef- 
fective impurity potential" is enlarged due to electron- 
electron correlation, which can be described by the GV 1 - 
method in Ref. [31]. For a reliable study of n in under- 
doped systems, it will be necessary to take account of 
residual disorders using the Gt^ 7 -mcthod. 

In summary, we studied thermal conductivity k in 
HTSCs. In the hole-doped shows a prominent 

"coherence peak" below T c , whereas it is absent in the 
electron-doped case. Based on the FLEX approxima- 
tion, such a contrasting behavior of k is well explained, 
although both YBCO and NCCO are pure d x 2. y 2-wa,ve 
superconductors. We do not have to assume a full-gap 
state in NCCO (such as d + is) to explain the absence of 
a coherence peak, which originates from the fact that the 
cold spot (line) in the normal state [~ (tt, 0)] is not on 
the nodal point (line) of the SC gap. The present study 
will open the way for the theoretical study of k in various 
interesting unconventional superconductors. 
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